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We here discuss a new model of 7-ray bursts (GRBs) based on differentially rotating strange 
stars. Strange stars in this model and differentially rotating neutron stars in the Kluzniak- 
Ruderman model can produce extremely relativistic, variable fireballs required by GRBs and 
then become millisecond pulsars. The effect of such pulsars on expansion of the postburst 
fireballs through magnetic dipole radiation is studied. We show that these two models can 
explain naturally not only various features of GRBs but also light curves of afterglows. 

PACS numbers: 97.70.Rz, 12.38.Mh, 26.60,+c, 97.60.Jd 

7-ray bursts (GRBs) has been one of the greatest mysteries in high-energy astrophysics 
for about thirty years. The observations of GRBs and a few afterglows [1-3] show that 
their sources are at cosmological distances. These observations require that any energy- 
source model should satisfy the following general features: (i) The model should produce 
an extremely relativistic fireball containing an amount of energy rs 10 51 -10 52 ergs [4]. 
(ii) The high Lorentz factor of the fireball implied by the short variability and the non- 
thermal spectra [5] requires that the fraction of contaminating baryons be less than 1% of 
the explosive energy (viz., < 1O _5 M ). (iii) The fireball is highly variable (required by the 
multi-peak feature of the light curves of GRBs) and it should last for the GRB duration 
(typically a few dozen seconds) , implying that the source cannot release energy only once 
[4,6]. (iv) The fireball should be a rare event occurring about once per million years in 
a galaxy [4]. (v) Finally, the model should account for the decline- rise-decline behavior 
of the R-band light curve of the afterglow from GRB 970508 [7,8]. In this Letter, we 
combine features of some existing pulsar GRB models, build upon them, and synthesize 
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a new model of GRBs based on differentially rotating strange stars. 

Recently Kluzniak and Ruderman [9] argued differentially rotating neutron stars as an 
origin of GRBs: In a differentially rotating neutron star, internal poloidal magnetic field 
(B) will be wound up into a toroidal configuration and linearly amplified as one part of 
the star rotates about the other part. Only when it increases up to a critical field value, 
Bf, will the toroidal field be sufficiently buoyant to overcome fully the stratification in 
neutron star composition. And then the buoyant magnetic torus will be able to float up 
to break through the stellar surface. Reconnection of the newborn surface magnetic field 
will lead to a quickly explosive event with a large amount of energy, which could be a 
sub-burst of a GRB. This model accounts well for the typical fluence in each observed 
sub-burst (peak), for the number of peaks, for the time interval between peaks, and for 
the rapid rise times and variability. This model may result in low-mass loading baryonic 
matter. 

Now let's discuss a new model of GRBs based on rotating strange stars in detail. 
Strange stars have been widely studied [10,11]. However, two objections against the exis- 
tence of isolated strange stars (with mass ~ 1.4M ) result from astrophysical arguments. 
On the one hand, it has been argued that the disruption of a single strange star can con- 
taminate the entire galaxy, and thus all neutron stars have in fact converted to strange 
stars [12]. Furthermore, if strange stars can be created directly in supernovae, then some 
strange stars must have more massive compact companions, and merging of such binaries 
will lead to the Madsen- Caldwell- Friedman effect [13]. This result conflicts with the post- 
glitch behavior of pulsars which has been described well by the neutron-superfluid vortex 
creep theory, rather than by the strange-star models [14] . On the other hand, an isolated 
neutron star with a stiff equation of state (EOS) and with mass of ~ 1.4M cannot con- 
vert to a strange star, because its central density is lower than the deconfinement density 
of neutron matter ~ (6-9)p (here p is the nuclear-matter density) [15]. 

However, there are two scenarios by which strange stars are born, and which don't 
conflict with the above two arguments. First, it has been proposed [16] that when neutron 
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stars in low-mass x-ray binaries accrete sufficient mass, they may convert to strange stars. 
The conversion should proceed in a detonation mode (at the speed of sound) [17], and 
thus the timescale for this process is about 0.1 ms. The direct criterion for the conversion 
is that the total masses of accreting neutron stars exceed ~ 1.8M Q if the EOS at high 
density is rather stiff. Fortunately, the theoretical analysis shows that the amount of 
matter accreted by the radio pulsars in many binary systems exceeds ~ 0.5M Q [18]. If 
this is true, some of these pulsars should have masses over 1.8M Q . The observation of 
kHz quasi-periodic oscillations in eight low-mass x-ray binaries and the analysis of high- 
resolution optical spectroscopy of Cyg X-2 indicate that the masses of neutron stars in 
these systems are larger than 1.8M [19]. 

We consider an alternative scenario. It is usually thought that due to gravitational 
radiation a neutron-star binary will eventually merge to a black hole. This results from 
the argument that the EOS for neutron matter at high density is soft. But, the EOSs 
in neutron stars are thought to be moderately stiff to stiff for the following reasons: (i) 
Soft EOSs at high densities are ruled out by the postglitch recovery in four pulsars [20]. 
(ii) The high masses (> 1.8M ) of the neutron stars in a few low-mass x-ray binaries 
also rule out soft EOSs [19]. (iii) The soft EOS such as kaon condensation doesn't seem 
to occur in stable neutron stars [21]. If the EOS for neutron matter is sufficiently stiff, 
therefore, the post-merger objects of Hulse- Taylor-like binaries may be massive neutron 
stars rather than black holes. The same outcome would be achieved if the initial masses 
of the merging neutron stars were low, e.g. M ~ 1M . According to the first scenario, 
these massive neutron stars will subsequently convert to strange stars. 

If strange stars form in some x-ray binaries by these scenarios, we now show that 
such stars are very difficult to be disrupted. The total number of strange stars in a 
galaxy is about 10 4 and their number in globular clusters (GCs) is N ss ~ 10 3 . The 
neutron star number density in GCs is n ns ~ 3 • 10" 59 (-R 9 /10kpc)~ 3 £ cm~ 3 , where R g 
is the galactic radius and £ ~ 10, and we have assumed that the total neutron-star 
number is ~ 10 s . The cross-section for a neutron star capturing a strange star is usually 
approximated by a ~ 4nG 2 M 2 /v 4 , where M is the mass of the neutron star and v its 
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velocity. Therefore, we find the rate for neutron star-strange star capture events in GCs, 
3? ~ avn ns N ss ~ 10- 17 (M/M o ) 2 (t;/500kms- 1 )- 3 ( J R 9 /10kpc)~ 3 (^/10) yr" 1 , and the rate 
outside GCs should be smaller. The actual rate for the strange stars to be first captured 
and then disrupted should be much less than 3ft due to the requirement of massive neutron 
stars whose number is rather small. 

For the above two scenarios, the pre- conversion neutron stars all have periods of the 
order of 1 ms because of accretion-induced spin-up in low-mass x-ray binaries and syn- 
chrotronization during the coalescence of neutron-star binaries. Since the moments of 
inertia always decrease during the conversion, the resulting strange stars will rotate at 
periods < 1 ms. Furthermore, differential rotation may occur in the interiors of these 
newborn strange stars due to the fact that the density profile of a strange star is much 
different from that of a neutron star with the same mass [11]. According to the basic idea 
of the Kluzniak-Ruderman model, such differentially rotating strange stars could lead to 
a series of sub-bursts of GRBs. 

Our model can account for the following features. First, sufficient energy. The total 
available energy released in our model should be the differential rotation energy, viz., 
E = iT2 2 /2 = 5 • 10 51 ergs Iu^l, where / is the effective moment of inertia of differential 
rotation (Z44 = 7/10 44 gcm 2 ), and Q the effective difference in angular velocity of rotation 
(f2 4 = Q/lO A s~ 1 ). We here estimate the critical field based on B"j/(8n) = fpcj., where / is 
the dimensionless parameter accounting for the stratification (~ 1% for massive strange 
stars), p the typical mass density of the strange star, and c s the speed of sound (c s = c/ 
for strange matter). Thus, we obtain B f rs 2 • 10 17 (p/4 • 10 14 g cm~ 3 ) 1/2 G. 

Second, very-low-mass baryon contamination. Because the strange star has a thin 
crust (~ 1O~ 5 M ) [10,11], the resulting fireball in our model is almost free of baryons, 
alleviating the baryon-contamination problem. This point was raised by Haensel et al. 
[22], who considered merging of strange stars as an origin of GRBs. Of course, this 
problem can also be avoided in many models [23-25], where power is extracted from 
rapidly spinning neutron stars or black holes with strong magnetic fields by Poynting 
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flux, which are later converted to 7-rays and lower-frequency photons. But, it is unclear 
whether these models can also explain the spiky light curves of GRBs. 

Third, temporal property. Following [9], we estimate the time interval between sub- 
bursts, Tj = (27r / 'Q) (B f I B) lOs-B-j^fiJ 1 , and the total duration, r d m 60s-B 13 1 f2 4 , where 
B13 is in units of 10 13 G. The internal field could be generated by either field-freezing or 
by a dynamo. Duncan and Thompson [23] have shown that whether a dynamo operates 
due to convection driven by a large neutrino flux is dependent on the Rossby number 
Nr = P/r conv , where P is the stellar period and r conv the convection overturn timescale 
at the base of the convection zone. If Nr < 1, an efficient dynamo can result and the 
(poloidal) field grows exponentially to B ~ 10 15 G; if Nr > 1, a dynamo is suppressed. 
It means that if a dynamo has operated, the duration is r d xs 0.6 s(5i3/10 2 ) _1 fi 4 , leading 
to short bursts; but if not, long bursts will form. This provides an explanation for the 
bimodal duration distribution of GRBs. This distribution has also been proposed to be 
associated with submillisecond neutron stars through gravitational instability [25]. 

Fourth, low burst rate. It has been argued [16,26] that the rate of conversion of accreting 
neutron stars in low-mass x-ray binaries, ~ 10~ 6 /yr per galaxy, is consistent with the 
observed GRB rate. In addition, the rate of coalesence of Hulse- Taylor-like binaries has 
been shown to correspond closely to the observed GRB rate [27]. 

Both the Kluzniak-Ruderman model and our model can produce an extremely rela- 
tivistic, variable fireball, which will lead to internal shocks [28] and subsequent external 
shocks [29]. A GRB will be produced once the kinetic energy is dissipated and radiated 
as 7-rays through synchrotron or possibly inverse-Compton emission from the accelerated 
electrons in these shocks. At the same time, the strange star is becoming a pulsar. There- 
fore, two common by-products of these models are a strongly magnetic millisecond pulsar 
and a postburst relativistic fireball. It is natural to expect that the central pulsar affects 
the evolution of the postburst fireball and the afterglow from this fireball. Can this effect 
explain the observed afterglows? Yes. Here we would provide a case in which this effect 
can explain well the optical afterglow of GRB 970508. 
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At the center of the fireball, the pulsar loses its rotational energy through magnetic 
dipole radiation, whose power varies with time as L(t) oc (1 + t/T)~ 2 , where t is one 
measure of time in the burster's rest frame, and T the initial spin-down timescale. For 
t < T, L can be thought as a constant; but for t^>T,L decays as oc t~ 2 . 

The pulsar radiates electromagnetic waves with frequency of u = 2tt/P ~ 10 4 s _1 . 
These waves are absorbed by the shocked interstellar medium (ISM) because the plasma 
frequency of the shocked ISM is much higher than uj. This implies that the pulsar pumps 
continuously its rotational energy into the shocked ISM. We assume that the expansion 
of the fireball in uniform ISM is relativistic and adiabatic. At a time t, the shocked ISM 
energy is given by 47rr 2 (r/47 2 )7 2 e', where r « ct is the blastwave radius, 7 the Lorentz 
factor of the fireball, and e' = A^ 2 nm p c 2 the shocked ISM energy density in the comoving 
frame (where n is the electron number density of the unshocked ISM) [30]. This energy 
should be equal to the sum of half of the initial energy (E/2) and the energy which the 
fireball has obtained from the pulsar based on energy conservation: 



where f3 — (1 — I/7 2 ) 1 / 2 . Because the Lorentz factor of the fireball at the initial stage 



pulsar is estimated as 4"f 2 E/L, measured in the burster's rest frame. The corresponding 
observer-frame timescale (r) is equal to this timescale divided by 27 2 , viz., r = 2E/L. 



on the stellar magnetic field. Please note that Blackman and Yi [31] recently discussed 
the dissipation of the rotational energy of a pulsar by large amplitude electromagnetic 
waves, which are directly converted to an afterglow, but in our model, afterglows result 
from shocks absorbing such waves. 

We consider below the synchrotron radiation from the accelerated electrons behind the 
shock, and assume a power-law distribution of the electrons with index p. We analyze 
evolution of the afterglow: (i) At the initial stage, viz., observer-frame time t e <C r, the 
second term on the right side of equation (1) can be neglected and the Lorentz factor of the 




(1) 



decays as 7 oc t 3 / 2 , the timescale in which the fireball has obtained energy ~ E/2 from the 



We assume r <C T, or P < 1.6ms(£ , 5 i/4) 1 l 2 l]!?. Clearly, this inequality is independent 
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fireball decays as 7 oc t® 3 ^ 8 . In this case, we have found that the fireball will rapidly go into 
the slow cooling phase, and the radiation at optical band will come from those radiative 
electrons a few hours later after the initial burst. This conclusion is similar to that of 
the recent studies [32,33]. Therefore, the optical flux from the fireball F v oc v~ p / 2 t®~ 3p ^ 4 . 
The result p ~ 1.0 inferred from the observed spectrum index in two days following the 
initial burst [8] leads to the R-band flux oc t© 1 ^ 4 . This satisfactorily accounts for the slow 
decline of the R-band light curve in ~ 8 hours after the burst [8]. 

(ii) For T > t ffi 3> r, the term E/2 in equation (1) can be neglected and the Lorentz 
factor of the fireball decays as 7 oc t^ 1 ^ 4 . In this case, the comoving-frame equipartition 
magnetic field decreases as B' oc 7 oc t© 1 ^ 4 , and the synchrotron break frequency drops 
in time as u m oc ^ 3 B' oc t ffi 1 . At the same time, since the comoving electron number 
density is n' e oc 7 oc and the comoving width of the emission region Ar' ~ r/7 oc t 3 / 4 , 
according to Meszaros & Rees in [7], the comoving intensity l' v oc n' e B'&r' oc 4 /4 - So 
the observed peak flux density increases in time based on F Vm oc t%l 5 I' Um oc t e . For 
a power-law distribution, electrons with different Lorentz factors may have different ra- 
diative efficiencies. Sari et al. [32] have given a critical value, 7 C oc 7 _3 ^0 1 oc t© 1 ^ 4 , 
above which cooling by synchrotron radiation is significant. The corresponding syn- 
chrotron frequency decays as u c oc 7-B'7 2 oc t^ 1 . From [32], we find the optical flux 
F v = (i/ c /i/ m ) _(p - 1)/2 (i/ /iy c )-P/ 2 F Um oc z/- p / 2 4 2 " p)/2 , and therefore, p « 1.0 further leads to 
the result that the optical flux increases with time, being in agreement with the subsequent 
rise of the R-band light curve [8]. 

(iii) For i e ^> T, the power of the pulsar due to magnetic dipole radiation rapidly 
decreases as L oc t~ 2 , and the evolution of the fireball is hardly influenced by the stellar 
radiation. As in stage (i), the optical flux is F v oc v Because p ~ 2.2, inferred 
from the spectrum index observed two days later after the burst, the R-band flux oc £^ L2 , 
which is also well consistent with the observations [8] . 

The effect of millisecond pulsars on afterglows can also be shown schematically in Figure 
1. We further give constraints on stellar parameters. According to the definitions of r 
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and T, we obtain the stellar period, P = 0.6ms( J B 5 i/4)- 1 / 2 (T/2d)- 1 / 2 (r/0.4d) 1 / 2 / 4 1 5 /2 , and 
the magnetic field strength B 13 = 3.1(E 51 /4)- 1 / 2 (r/2d)- 1 (r/0.4d) 1 / 2 / 4 5^ 3 , where E S1 = 
S/10 51 ergs, I45 is the moment of inertia in 10 45 gcm 2 and Rq the stellar radius in 10 6 cm. 
The observations indicate r ~ 0.4 days and T ~ 2 days [8]. If taking (E^i/^^^I^Rq 3 ~ 
1, we have Bi 3 ~ 3.1, and the total duration of the main burst (to) is estimated to be about 
20 seconds. This is in excellent agreement with the observation by the BeppoSAX satellite 
[34]. In addition, we find P ~ 0.6 ms, implying that the pulsar has a submillisecond period. 
What we want to point out is that such a period for the formation scenario of a strange 
star in a low-mass x-ray binary may be due to preconversion accretion-induced spin-up 
and due to decrease of stellar moment of inertia during conversion, and a magnetic field 
~ 10 13 G may result from preconversion crustal plate motion and/or from a dynamo in 
the strange star. 

If a GRB results from magnetic reconnection on the surface of a differentially rotating 
strange star or neutron star, we cannot clearly determine the region in which such a 
process takes place. A beaming of the emission might form if this process occurs on a 
small area of the stellar surface as solar flares. This might provide a solution to the 
large-energy problem of GRB 971214. 

Finally, it is interesting to note that two other kinds of high-energy transients, the soft 
7-ray repeaters and the hard x-ray burster (GRO J1744-28), have also been described 
well by the strange-star models in which the stellar crust-breaking leads to bursts [35]. 

We thank the referees for their enlightening comments and valuable suggestions. This 
work was supported by the National Natural Science Foundation of China. 
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FIGURE CAPTION 

FIG. 1. Afterglow light curves with a pulsar (solid line) and without a pulsar (dashed 
line) . 
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